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We study how dephasing affects the distribution of the dc current pumped through a chaotic 
quantum dot. We introduce dephasing by the addition of a voltage probe to the quantum dot, 
treating both the case of controlled dephasing (when the voltage probe is coupled to the dot via a 
ballistic point contact with a conductance that is increased stepwise) and intrinsic dephasing (for 
which the voltage probe serves as a phenomelogical model). While dephasing eventually suppresses 
the dc current through the dot, we also find that, for a quantum dot with single-mode point contacts, 
a small amount of dephasing actually decreases the likelihood of a zero pumped current. 

PACS numbers: 72.10.Bg, 73.23.-b 



I. INTRODUCTION 

Quantum mechanical interference is a key player deter- 
mining sample-specific properties of small metal or semi- 
conductor particles. For closed systems, quantum inter- 
ference fixes the precise positions of energy levels and the 
microscopic details of wavefunctions; for open syste 
it governs, e.g., the fluctuations of the conductance 
In all these cases, interference provides a (sometimes 
large) correction added to a nonzero background. Such 
a background is absent in a so-called "quantum electroa 
pump" , which Jias. received considerable experimental^ 
and theoreticaloil3 attention recently. This requirement 
of phase coherence for the observation of a pumped cur- 
rent makes a quantum pump a sensitive instrument to 
study quantum interference in quantum dots. 

In the experimental proposal of Ref. [|, a quantum 
electron pump is made of a semiconductor quantum dot 
connected to two electron reservoirs. A dc current is 
generated by periodically varying two gate voltages that 
characterize shape of the dot. (Alternatively, one may 
pump electrons by periodic variation of the Fermi en- 
ergy, magnetic flux, etc.) The electron pump is called 
"quantum" , because the variation pertains to the quan- 
tum interference only, not to the classical dynamics of 
the dot: In a trajectory-based picture, the variations af- 
fect phases of the trajectories, not their weights. For 
this reason, this "quantum" electron pump is fundamen- 
tally different from electron pumps or turnstiles that rely 
on Coulomb., charging physics and do not require phase 
coherence .Eailil In the quantum pump, the precise value 
of the pumped current and its direction depend on micro- 
scopic details and vary from sample to sample. There- 
fore, one needs to address the probability distribution 
P(I) of pumped currents for an ensemble of quantum 
pumps. Theoretical work has focused on electron pumps 
in a chaotic quantum dot, for which this distribution can 
be calculated using random matrix theory.ffi3&[L3 For a 
pump with single-channel point contacts, it was found 



that P(I) has a cusp at / = and algebraic tails,Q ; 
for many-channel point contacts, P(I) is Gaussian.! 

In this paper, we study how the probability distribu- 
tion of the pumped current is affected if phase coherence 
in the quantum pump is gradually destroyed. As in Rcfs. 
p|,p|~^0|, we consider the case of a quantum pump that 
consists of a chaotic quantum dot. A controlled way to 
destro¥-phase coherence is to couple the dot to a voltage 
probe.E-3 While not drawing a net current, the voltage 
probe absorbs and reinjects electrons without any phase- 
relationship, thereby destroying phase coherence. The 
dephasing rate 74, can be tuned by varying the conduc- 
tance of the point contact connecting the quantum dot 
and the voltage probe. While a voltage probe can be used 
as an external source to controllably increase the dephas- 
ing rate in the quantum pump, it can also be used as a 
model for intrinsic dephasing processes in the pump itself, 
such as electroa electron interactions or two-level systems 
in the dot.Ej EO A similar approach was taken in a recent 
paper by Moskalets and Buttiker, who studied the ef- 
fect of a voltage probe on the current pumped through 
an electron pump consisting of a one-dimensional wice 
with two tunable tunnel barriers of oscillating strength.E£l 
However, unlike the quantum dot geometry studied here, 
the pump of Ref. [l8] can also operate as an electron pump 
in the absence of phase coherence. The variance of the 
pumped current in the case of many-channel point con- 
tacts was previously obtained by Shutenko et al. using a 
different model B 

In Sec. we use the scattering approach to derive a 
formula for the pump ed current in the presence of the 
voltage probe. In Sec. Ill we then study the distribution 
of the pumped current for the cases of controlled dephas- 
ing and intrinsic dephasing. We focus on the cases of a 
quantum dot with single channel point contacts and with 
multi-channel point contacts. We conclude in Sec. [V. 
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II. PUMPED CURRENT IN THE PRESENCE OF 
DEPHASING 

A schematic picture of the system is shown in Fig j^. 
It consists of a quantum dot connected by ballistic point 
contacts to three electron reservoirs, labeled 1, 2, and 
3. The point contacts allow N\, N 2 and A3 propagating 
channels at the Fermi level, respectively. The total num- 
ber of channels in all point contacts will be denoted as 
N = Ni + N 2 + A 3 . The reservoirs 1 and 2 are held at 
the same voltage V\ = V 2 = 0, while reservoir 3 serves 
as a voltage probe: the time-dependent voltage Vs(t) is 
adjusted such that the current Is(t) = at all times. 

Two external parameters Xi(t) and X2(t) that deter- 
mine the shape (or other characteristics) of the quantum 
dot are varied periodically with frequency w. As a result 
of this periodic variation, a dc current 1 = 1% = —I2 will 
flow between reservoirs 1 and 2. We will now evaluate 
this pumped current in the presence of the voltage probe. 
Starting point of our evaluation is the relatiorfl between 
the pumped current and the emissivity edn(m)/dXj, the 
charge that leaves the dot through contact m = 1,2 as 
the parameter Xj is varied by an amount dXj, 
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Here A is the area in (Xi, A2)-space enclosed by the pa- 
rameters [Xi(t) , X2(t)] in one cycle. Equation (p is valid 
for the adiabatic regime, ujTd <C 1, where ~ h/NA is 
the dwell time of the quantum dot, A being the mean 
single-particle level spacing. 

In the absence of the voltage probe, and for ballistic 
point contacts between the dot and reservoirs 1 and 2, 
so that (effects of charge quantization in the dot can be 
ignored,Ej the emissivity edn(m) /dXj ia-j-elated to the 
A x A scattering matrix S of the dot asc3 
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With voltage probe, charge can either leave the quantum 
dot directly through the contacts 1 or 2, or via inelastic 
scattering in reservoir 3. Hence, we can write the emis- 
sivity dn(m)/dX as a sum of an elastic and an inelastic 
contribution, 
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where the elastic contribution (dn(m) / dX) c \ is still given 
by Eq. (J2J) and the inelastic contribution reads 
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is an element of the 3x3 condutance matrix of the dot. 
Substitution of Eqs. (§)-(§) into Eq. (0) then yields the 
current formula 
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The first contribution idir represents the charge that ex- 
its the dot either directly or indirectly via the reservoir, 
while the second contribution / rcct is an additional con- 
tribution to the dc .current that arises from rectification 
of the voltage p3(t).E3 One verifies that in the limit where 
the third reservoir is decoupled from the quantum dot, 
-frcct vanishes, while the expression for Tdir approaches 
that for current in the absence of dephasing (see Ref. ||). 

In the next section we proceed to calculate the sta- 
tistical distribution of the current I for an ensemble 
of chaotic quantum dots. We restrict ourselves to the 
case of small pumping amplitudes X\ — 5X\ sm.{iot), 
X 2 = 5X 2 sin(wi + <p), for which the current / is bilinear 
in 8X\ and 5X 2 , 



I = — eui sin <pSXi5X 2 , 
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We will consider both the case of "controlled dephas- 
ing" , corresponding to a real voltage probe coupled to 
the quantum dot via a ballistic point contact, and that 
of intrinsic dephasing, which is modeled by a voltage 
probe coupled to the quantum dot via a wide tunneling 
contact .Ell 

In the absence of dephasing, the pumped current is not 
symmetric undet reversal of a magnetic field through the 
quantum dot .BO One easily verifies that in the presence 
of dephasing the current / remains asymmetric under 
reversal of the magnetic field, even in the presence of 
strong dephasing (A 3 large). 



III. DISTRIBUTION OF THE PUMPED 
CURRENT IN THE PRESENCE OF DEPHASING 

For an ensemble of chaotic quantum dots the distri- 
bution of the scattering matrix and its derivatives is 
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FIG. 1. A quantum dot with two parameters X\ and X2 
which describe a deformation of its shape. As X± and X2 are 
varied periodically a dc current I — I\ — —I2 flows from the 
right to the left reservoir. The lead on the bottom, through 
which no current flows (^3 = 0), has a variable number of 
channels A3. 

known.0 It is most conveniently expressed through the 
parameterization 



S = UU', 
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where U and V are N x N unitary matrices and Qj is an 
N x N hermitian matrix. In the presence of time-reversal 
symmetry, U' = U T and Qj = Qj. We also introduce 
the matrix Q e that parameterizes the energy- derivative 
dS/de, 

OS 2ir 

For ballistic point contacts to reservoirs 1, 2, and 3, the 
matrices U and V are uniformly distributed in the uni- 
tary group, independently of Q±, Q2, and Q £ , while the 
joint distribution of the matrices Qx, Q2, and Q e is given 
by 



Poc(detQ e )- N / 2 -W N+2 -Ve(Q) 
x exp 
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where <d(Q) = 1 if all eigenvalues of Q are positive and 
O(Q) = otherwise and = 1 (2) in the presence (ab- 
sence) of time-reversal symmetry. Note that at fixed Q e 
the distribution of Q\ and Q2 is Gaussian with a width 
set by Q £ . 

We will now study the distribution of the pumped dc 
current for controlled dephasing and instrinsic dephasing. 
In each case we will first study the variance and then the 
full distribution. 



A. Controlled Dephasing 

The effect of dephasing on the distribution of the 
pumped current can be studied in a controlled setting 
by increasing the number A^ of open channels in the 
voltage probe one by one. The third reservoir serves as 
a true voltage probe if the charge relaxation rate of that 
reservoir is much larger than the frequency u>, so that 
the voltage V^{t) can adjust essentially instantaneously 
to balance any current Ia(t) flowing into or out of that 
reservoir as a result of the pumping action on the quan- 
tum dot. 

A voltage probe connected to the dot via a ballistic 
point contact with N3 channels gives rise to a dephasing 
rate 
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We now fix N\ and N2 and calculate the distribution of 
the pumped current as a function of A%. We consider 
the cases Ai — N2 = 1 of single channel current-carrying 
leads and N% , N2 3> 1 of many-channel current-carrying 
leads. In the former case, we have calculated the first two 
moments of the distribution P(i) of the dimcnsionlcss 
current i analytically, using the technique of Ref. [2^ to 
perform the integrations over the matrices U and U'. 
The average current (i) was thus found to be zero, while 
the variance (i 2 ) decreases with increasing N3; 
divergent for A^ = for (3 = 2 and for A^ < 6 for 
The results of this calculation are shown in Fig. 
large A3, we find 
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Details of the calculation can be found in Appendix [A|. 
[Equation ( |Io| ) is valid for arbitrary Ni and A2, as long 
as A3 ^> 1.] For A^ ^> 1, the contribution of the recti- 
fication current i rcc t to (i 2 ) is proportional to A^ 4 and 
hence negligible. (A similar conclusion for the rectifica- 
tion contribution was reached in Ref. [l| for the compar- 
ison of measurements of pumped current and pumped 



voltage.) As shown in Fig. |2|, in general, 



is larger 



with time-reversal symmetry than without, though the 
difference between the two variances vanishes for large 
dephasing rates. 

We obtained the full distribution P(i) for N\ = N2 = 1 
using Monte-Carlo integration with the distribution (jg). 
The result is shown in Fig. ||for various values of A3. De- 
tails of the numerical procedure are outlined in App. 
Note that the distributions are symmetric around i = 0. 
In the absence of dephasing (A3 = 0), P(i) has a cusp 
at i = 0, which is smoothed out for A3 > 0. The distri- 
bution P(i) has algebraic tails with powers that increase 
with increasing dephasing rate. In the presence of time- 
reversal symmetry the tails fall off with a power smaller 
than 3 as long as ^3 < 6, corresponding to a diver- 
gent second moment (i 2 ) (cf. Fig. 0). Interestingly, while 
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FIG. 2. Variance of the pumped current with one open 
channel in each of the current carrying leads (iVi = N2 = 1) 
as a function of the number of open channels in the third lead 
(N3). The presence (absence) of time-reversal symmetry is 
denoted by open circles (closed circles). The dashed line is 
the variance in the asymptotic limit Ni = N% = I, N3 — > 00. 
The variance in the presence of time- reversal symmetry is 
infinite for N3 < 6. 



dephasing decreases the second moment (i 2 ) for single- 
channel point contacts iVi = N2 = 1, a small amount 
of dephasing also decreases the probability P(i = 0) of 
finding no pumped current at all. This can be under- 
stood physically, since the cusp in the distribution P(i) 
at zero pumped current arises from a the occurence of 
destructive interference, which is as much suppressed by 
dephasing as the constructive interference responsible for 
large values of % and the tails of the distribution. For 
JV3 ^> 1, the distribution P(i) approaches a Gaussian. 

For multichannel point contacts in the current-carrying 
leads, N\,N2 3> 1, the distribution P(i) is a Gaussian 
centered around i = 0. To calculate the second mo- 
ment (i 2 ), we parameterize the scattering matrix S and 
its derivatives as in Eq. (]?]), perform the Gaussian inte- 
gration over Q\ and Q2 with the distribution (||) at fixed 
Q e , integrate over the unitary matrices U and V using 
the diagrammatic technique of Ref. ^2] and finally inte- 
grate over the eigenvalues Tj, j = 1, . . . , N, of Q £ , see 
appendix [A| for details. The result is given by Eq. dl" 
above. 
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FIG. 3. The probability distribution of the pumped cur- 
rent without (a) and with (b) time-reversal symmetry, for 
the case Ni — N2 = 1 of single-channel point contacts in 
the current-carrying leads and as a function of the number 
of channels A^3 in the voltage probe. The dashed line shows 
the distribution without dephasing (JV3 = 0). The solid lines 
show the distribution with N3 = 1, (lowest curve at i = 0) 
2, ... , N 3 = 10 (highest curve). 



B. Intrinsic dephasing 



The voltage probe can also be used as a model for in- 
trinsic dephasing in the quantum dot. Sources of dephas- 
ing may be, e.g., electron-electron interactions or interac- 
tions with an external bath of photons and/or phonons. 
In such a case the source of decoherence is delocalized 
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throughout the dot. This situation is well modeled by 
a voltage probe with a tunnel barrier with transmission 
probability r 3 « 1 and many channels A 3 3> 1, such 
that the product 
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is kept fixed. 

In order to find the distribution of the pumped current 
in this case, we need the distribution of the scattering ma- 
trix S and its derivatives when the contact to the third 
reservoir contains a tunnel barrier. This problem can be 
solved by a statistical mapping which connects the scat- 
tering matrix S to a scattering matrix So that is taken 
from an ensemble appconciate for a quantum dot with 
ballistic point contacts,E3a 



s = x/T - r- Vf 
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for 



where T is a diagonal N x N matrix with Tj 
index j corresponding to the current-carrying leads 1 and 
2 and Tjj = r 3 for index j corresponding to the voltage 
probe. The distribution of So and its derivatives is as 
described in Eqs. (0)~(||) at the beginning of this section. 
We can then find the distribution of S from Eq. (|lj) and 
the distribution of its derivatives upon differentiating Eq. 

We have calculated the variance (i 2 ) and the full distri- 
bution of the current for Ni = N2 = 1 using Monte-Carlo 
integration with the above distribution. The results are 
shown in Figs. | and ||. We note that intrinsic dephasing 
cuts off the tails of the distribution P(i). This is in con- 
trast to the case of controlled dephasing by a few-channel 
voltage probe with ballistic point contacts, which merely 
replaces the algebraic tail of the distribution at zero de- 
phasing by another, faster decaying, algebraic tail. As 
for controlled dephasing, the variance of the pumped cur- 
rent is larger with time-reversal symmetry than without, 
though for intrinsic dephasing the difference is smaller 
than in the case of controlled dephasing (cf. Fig. |^). At 
large dephasing rates, the dependence of (i 2 ) on the pres- 
ence or absence of time-reversal symmetry vanishes. Also 
note that the probability to find zero pumped current 
initially decreases when the dephasing rate is increased, 
although the effect is not as strong as in the case of con- 
trolled dephasing. For large dephasing rates, both intrin- 
sic and controlled dephasing yield the same distribution 
P(i). 

For large Ni and N2 the distribution P(i) is again 
Gaussian, with zero mean and with variance given by 
Eq. (10) above with A3 replaced by 7. This result agrees 
with what was previously obtained by Shutcnko et al. 
using a different method n 
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FIG. 4. The variance of the pumped current in the case 
of intrinsic dephasing as a function of dimensionless dephas- 
ing rate 7 in the presence (dashed curve) and absence (solid 
curve) of time-reversal symmetry. The asymptotic result (dot- 
ted curve) is also shown. There is one open channel in each 
lead (Ni =N 2 = l). 
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FIG. 5. The probability distribution of the pumped cur- 
rent with Ni = -/V2 = 1 in the presence of time-reversal sym- 
metry and intrinsic dephasing. The dashed curve is the proba- 
bility distribution in the absence of dephasing. The two solid 
curves have dimensionless dephasing rates 7 = .1 (highest) 
and 7=1 (lowest). The dot-dashed curves have dephas- 
ing rates 7 = 8 (lowest) and 7 = 15 (highest). The inset 
shows the probability distribution of the pumped current with 
Ni — N2 — 1 in the absence of time-reversal symmetry with 
the same parameters. The distribution with 7 = .1 is indis- 
tinguishable from the case without dephasing (7 = 0). 
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IV. CONCLUSION 

In summary, we have derived the current distribu- 
tion of electrons pumped adiabatically through a chaotic 
quantum dot, in the presence of a voltage probe. The 
voltage can either serve as a controlled source of dephas- 
ing, or as an effective description of intrinsic dephasing 
processes inside the quantum dot. For the case of a quan- 
tum dot with two single-mode current-carrying point 
contacts, the conductance distribution is non-Gaussian 
with algebraic tails. Remarkably, while dephasing shifts 
the weight of the probability distribution P(I) of the 
pumped current towards zero current, a small amount of 
dephasing actually reduces the probability P(I = 0) to 
find no pumped current at all. This can be understood if 
the probability to find zero pumped current is enhanced 
by destructive interference — note the cusp in P(I) at 
1 = — , which is then suppressed by dephasing. For a 
quantum dot with many-channel point contacts (as well 
as for a quantum dot with single-channel point contacts 
and a large dephasing rate), the current distribution is 
Gaussian. The width of the distribution decreases mono- 
tonically with increasing dephasing rate. For dephasing 
rates 7^ much larger than the escape rate 7 ~ (Ni+N 2 )A 
to the reservoirs, the r.m.s. current decays oc (7/7^) 3 ^ 2 
(see also Ref. ||). 

Not only dephasing, but also thermal smearing can 
reduce the size of the pumped current. Thermal smear- 
ing has been considered by Shutenko et aZ.,B who found 
that, without dephasing, the r.m.s. current decays oc 
(7/T) -1 / 2 for temperature T 3> 7. Experimentally, the 
dephasing rate 7^ has been found to increase oc T (or 
faster), which implies that, for such a dephasing rate, 
dephasing is the more effective mechanism to reduce the 
current pumped in a quantum pump. 

It is interesting to compare the quantum pump to 
a (quantum) rectifier. The latter device generates a 
dc current in response to an a.c. bias voltage. For 
a setup similar to the quantum dot pump considered 
here, the variance of-, the rectified current is propor- 
tional to (dG/dX) 2 E3 For^-large dephasing rates, typi- 
cally (dG/dX) 2 ~ (7/70) 3 ,E3 so that the r.m.s. rectified 
current decays oc (j/jip) 3 ^ 2 - As in an experimental real- 
ization a true quantum pump may coexist with a recti- 
fier (the a.c. voltages arising frorrudisplacement currents 
and parasitic capacitive coupling) ,EZI we conclude that de- 
phasing does not change the relative importance of one 
mechanism over the other. 
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APPENDIX A: 

In this appendix we outline the calculation of the mo- 
ments (i 2 ) in the presence of controlled dephasing. A 
detailed report of the calculations will appear in Ref. p^. 

In order to calculate the average (i 2 ), we need to per- 
form an average over the N x N matrices U, U', Q%, and 
Q2 defined in Eq. (^). The matrices U and U' are uni- 
formly distributed in the unitary group [LP = U T in the 
presence of time- reversal symmetry (TRS)], while the dis- 
tribution of the matrices Q\ and Q2 is given by Eq. (jg). 
First we perform the average over Q\ and Q2. Hereto, 
we parameterize Qi and Q2 as 



Qi = 



1,2, 



(Al) 



where Hi and H2 are hermitian N x N matrices (real 
symmetric in the presence of time-reversal symmetry) 
and $ is a complex (real) N x N matrix such thatEa 
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Substitution of Eqs. ( |Al| ) and (A2) into the distribution 
(0) shows that the elements of Hi and H 2 have a Gaus- 
sian distribution with zero mean and with variance 
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We first perform the Gaussian average over the ma- 
trices Hi and i?2- The resulting expression contains the 
matrix ^ in the combination Q e = ^i^ 1 ^^ 1 only. We 
decompose Q e in terms of its eigenvalue matrix f and a 
random unitary (orthogonal) matrix V, 



= VtVK 



(A4) 



and average over V using the diagrammatic method of 
Ref. (or its straightforward generalization to orthogo- 
nal matrices V in the case when time-reversal symmetry 
is present). The result of the average is an expression 
that depends on the scattering matrix S and the matrix 
of eigenvalues f only. 

In the absence of time-reversal symmetry, the result of 
this calculation is 
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while in the presence of time-reversal symmetry, we find 
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The contribution from i rect is a factor N smaller than 
idii for large N. Now, the large- N result, Eq. (|Io|), fol- 
lows easily employing the known density of dimcnsionless 
delay times, 
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where r± = {3± VE)/N, to integrate over f. 

For the case Ni = N 2 = 1, we were able to obtain 
exact expressions for («dir*rcct) and (i^ Bct ) for arbitrary 
N 3 and time-reversal symmetry is broken. These expres- 
sions, which contain up to a product of four traces involv- 
ing the matrix f, were too lengthy to be reported here, 
and will be published in Ref. The results for (i 2 ) are 
shown in Fig. ^. We summarize the main steps of the 
calculation below. Starting from Eq (A5) and similar 
expressions for the variance of i rcc t£j we integrate over 
the eigenvalues Tj of the dimensionless time-delay ma- 
trix Q £ and the scattering matrix S. To average over the 
Tj, j = 1, . . . , N, we introduce the dimensionless escape 
rates x n = l/r„ which are distributed according to the 
generalized Laguerre ensemblejEl 

P(xi, . . . ,x N ) (x Y\ (x n -x m ) 2 Y[ x^exp(x n ). (A8) 

n < rn n 

Since the r„ appear only in products of up to four 
traces, we need the marginal n-point distributions 
Rn(xi, . . . , x n ) for n = 1, . . . , 4 only, where 



Rn {x \ , ■ ■ ■ 5 X n *j 



x / dx n+1 ...dx N P(x 1 ,...,x N ). 
Jo 



We can find exact expressions of i?i , . . . , R4 using the as- 
sociated Laguerre polynomials (x) which are orthog- 
onal with respect to the weight function x N e~ x . in par- 
ticular, for the case of broken time-reversal symmetry, 



Rn(xx, ...,x n ) = det [K(x i ,x j )) i 



where 



K(xi,Xi) = Y^L^{x l )L^{x ] )e^ +x ^/ 2 (x iXj ) 



N/2 



1 = 1 



The average over 5* can be done using the polar decom- 
position of 5, 



S = 



\ f Vl - tU itf 



v 



u' 

it VT^tH J V v ' 1 



where u and u' (v and v') are 2x2 (N3 x N3) unitary 
matrices and t is an N3 x 2 matrix with all elements equal 



to zero except t r , 



it: 



1,2. In the presence 



of time- reversal symmetry, vl = u T and v' = v T . The 
parameters T± and T2 govern the escape rate into the 
voltage probe. The uniform distribution of U and U' in 
the unitary group yields the integration measure 



dS = m 



T 2 \ 2 (T 1 T 2 f 3 



' dudu' ' dvdv 1 ' dT\dT 2 . 



The average over S then reduces to an integral over T\ , T 2 
and one angle <f> uniformly distributed in the interval < 
(f> < 7r/2. In terms of these variables we can write 

G13 = Ticos 2 + T 2 sin 2 

G 23 = T lS in 2 + T 2 cos 2 
tvS S3 St 3 S S3 St 3 = N - 4 + (1 - T1) 2 + (1 - T 2 f 
XxS\ 3 S X3 S\ 3 S 33 = (1 - T^cos 2 ^ + (1 - r 2 )r 2 sin 2 0, 

where Sij denotes the ij block in the scattering matrix-S. 
(The last two terms appear in the expression for z rcct .Ej) 



APPENDIX B: 

In order to obtain the full distribution of the pumped 
current j, we numerically generated matrices U, U', Q\ 1 
and Q 2 according to the appropriate distributions. The 
matrices U and U' are uniformly distributed in the uni- 
tary group (U 1 — U T in the presence of time-reversal 
symmatuy). Using a parameterization in terms of Euler 
angles,Lj their generation is relatively straightforward. 
The numerical generation of matrices Q\ and Q 2 accord- 
ing to the distribution (|J) makes use of a trick that was 
inspired by Ref. which would like to explain below. 

We parameterize the N x N matrices Q±, Q 2l and Q £ 

as 



Qe 



1,2 



(Bl) 



where C is a complex N x 2N matrix [real N x (2N + 1) 
matrix in the presence of time-reversal symmetry] and 
Hi (i = 1,2) is a hermitian 2N x 2N matrix [real sym- 
metric (2N + 1) x (2N + 1) matrix in the presence of 



time-reversal symmetry]. In Eq. (Bl), the inverse C 1 is 



defined as the right-inverse, C*" 1 = C T (CC T )^ 1 . In Ref. 
|3j] it was shown that the matrix Q e has the distribution 
(g) if the elements of C are all chosen from a Gaussian dis- 
tribution with unit var ian ce. From there, one can show 
by substitution of Eq. (Bl) into Eq. (||), that the param- 
eterization (Bl) reproduces the correct distribution (^) 
for all three matrices Q e , Qi, and Q 2 if the elements of 
the matrices H\ and H 2 are all chosen f rom a Gaussian 
distribution with variance given by Eq. (A3) above. 
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